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Abstract 

Existence and uniqueness of global in time measure solution for a one dimensional non- 
linear aggregation equation is considered. Such a system can be written as a conservation 
law with a velocity field computed through a selfconsistant interaction potential. Blow 
up of regular solutions is now well established for such system. In Carrillo et al. (Duke 
Math J (2011)) |14) . a theory of existence and uniqueness based on the geometric approach 
of gradient flows on Wasserstein space has been developped. We propose in this work 
to establish the link between this approach and duality solutions. This latter concept of 
solutions allows in particular to define a flow associated to the velocity field. Then an 
existence and uniqueness theory for duality solutions is developped in the spirit of James 
and Vauchelet (NoDEA (2013)) |21| . However, since duality solutions are only known in 
one dimension, we restrict our study to the one dimensional case. 
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1 Introduction 



Aggregation phenomena of a population of particles interacting under a continuous interaction 
potential are modelled by a nonlocal nonlinear conservation equation. Letting p denote the 
density of cells, the so-called aggregation equation in space dimension writes 

5iP + V. ■ (a(V.iy * p)p) = 0, t>0, xgM^, (1.1) 

and is complemented with the inital condition p(0,x) = p*"*. Here W : — M is the 
interaction potential, and a : ^ is a smooth given function which depends on the 
actual model under consideration. 

This equation is involved in many applications in physics and biology. In the framework of 
granular media [21 [IS 125], a is the identity function, and interaction potentials are in the form 
W{x) = —\x\°' with a > 1. In plasma physics, the context is the high field limit of a kinetic 
equation describing the dynamics of electrically charged Brownian particles interacting with a 
thermal bath. This leads to consider potentials in the form W{x) = — \x\, and a = id as well (see 
e.g. |27|). Also, continuum mathematical models have been widely proposed to model collective 
behaviour of individuals. Then the potential W is typically the fundamental solution of some 
elliptic equation, and a depends on the microscopic behaviour of the individuals. We refer for 
instance to the well-known Patlack-Keller-Segel model where aggregation of cells is described 
by a macroscopic non-local interaction equation with linear diffusion |30] . More precisely, 
the swarming of cells can be described by aggregation equations where the typical interaction 
potential is the attractive Morse potential W{x) = |e~l^l [T5| EHl EH]- Such potentials also 
appear when considering the hydrodynamic limit of kinetic model describing chemotaxis of 
bacteria [HI [HI [20]. 

All the potentials mentioned above have a singularity at the origin, they fall into the context 
of "pointy potentials" (see a precise definition below), and it is well-known that in that case 
concentration phenomena induce the blow-up in finite time of weak solutions [ll[5l[6]. Thus 
the notion of solution breaks down at the blow-up time and weak measure-valued solutions for 
the aggregation equation have to be considered. Carrillo et al. |14) have studied the multidimen- 
sional aggregation equation when a = id in the framework of gradient flow solutions. Namely, 
equation (11. ip is interpreted as a differential equation in time, the right-hand side being the 
gradient of some interaction energy defined through the potential W. This idea, known as the 
Otto calculus (see fIE\ [55]). requires the choice of a convenient space of probability measures 
endowed with a Riemannian structure. Then, following [2], gradient flow solutions are inter- 
preted as curves of maximal slopes in this space. The authors obtain existence and uniqueness 
of weak solutions for ( II. ip in M^, A^ > 1 when a = id, the main problem being now to connect 
these solutions to distributional solutions. 

An alternative notion of weak solutions has been obtained by completely different means 
in the framework of positive chemotaxis in [2T]. Here equation (II. ip with W{x) = |e~'^'' can 
be obtained thanks to some hydrodynamic limit of the kinetic Othmer-Dunbar-Alt system |16) . 
The key idea is to use the notion of duality solutions, introduced in [9] for linear conservation 
equations with discontinuous velocities, where measure- valued solutions arise also. In that case, 
this allows to give a convenient meaning to the product of the velocity by the density, so that 
existence and uniqueness can be proved. When applying this strategy to the nonlinear case, it 
turns out that uniqueness is not ensured, unless the nonlinear product is given a very precise 
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signification, see for instance |10[ [T9] . In tlie case of cliemotaxis, it is provided by tlie limit 
of the flux in the kinetic model. Once this is done, existence and uniqueness can be obtained. 
An important consequence of this approach is that it allows to define a flow associated to 
this system. Then the dynamics of the aggregates (i.e. combinations of Dirac masses) can be 
established, giving rise to an implementation of a particle method and numerical simulations 
of the dynamics of cells density after blow-up (see also [22] for a numerical approach using a 
discretization on a fixed grid). The principal drawback of this method is its limitation to the 
one-dimensional case, mainly because duality solutions are not properly defined yet in higher 
space dimension. Thus the theory developped in [T4| is, up to our knowledge, the only one 
allowing to get existence of global in time weak measure solution for (11. ip in dimension higher 
than 2. Another possibility could be using the notion of Filippov flow |18| . together with the 
stability results in [7], to obtain a convenient notion of solution to (II. ip . thus following |31| . 

This work is devoted to the study of the links between these two notions of weak solutions 
for equation (II. ip . in the one-dimensional setting. As we shall see there is no equivalence strictly 
speaking for a general potential and a nonlinear function a. More precisely we first consider 
the same situation as in jl4j . that is a pointy potential W , and a = id. We adapt the proof of 
|21| to define duality solutions in this context, and choose a convenient space of measures to 
be compatible with the gradient flows. Then we prove that duality solutions and gradient flow 
solutions are identical (Theorem 14.11 below), thus answering the questions raised by Remark 
2.16 of [H]. 

Next, we investigate the nonlinear case, that is a 7^ id. Notice that additional monotonicity 
properties are required to ensure the attractivity of the dynamics. The results of [H] cannot 
be applied as they stand, the key problem is to define a new energy for which weak solutions of 
(II. ip are gradient flows. However, we are able to find such an energy only in the particular case 
W = —^\x\. On the contrary existence of duality solutions for (II. ip with a nonlinear function a 
can be obtained for more general potentials W, even if we cannot reach the complete generality 
of [H]. 

The outline of this paper is as follows. In the next Section, we introduce notations and 
recall the main results obtained in [T3] in the case a = id. A sketch of their proof is proposed. 
Section [3] is devoted to the duality solutions, and starts by recalling their original definition and 
main properties. Next, we turn to the nonlinear setting, and define precisely the velocities and 
fluxes that allow to state the existence and uniqueness results both for a = id and a 7^ id. The 
case a = id is treated in Section HI existence and uniqueness for duality solutions are proved, 
together with equivalence between gradient flows and duality solutions. Finally in Section [5] we 
investigate the case a 7^ id, where general equivalence results no longer hold. 

2 Gradient flow solutions 
2.1 Notations and definitions 

Let Co(Y, Z) be the set of continuous functions from Y to Z that vanish at infinity and Cc{Y, Z) 
the set of continuous functions with compact support from Y to Z. All along the paper, we 
denote MiocC^^) the space of local Borel measures on R^. For p G MioM^) we denote by 
IpKR-'^) its total variation. We will denote A^fe(M^) the space of measures in AiiodR'^) with 
finite total variation. From now on, the space of measure-valued function A^f,(M^) is always 
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endowed with the weak topology a{Mb, Co). We denote Sm ■= C([0, T]; Mb(^^) -a{Mb, Co)). 
We recall that if a sequence of measure (yUn)nGN in -Mbi^^) satisfies sup^gf^ |/i„|(]R^) < +00, 
then we can extract a subsequence that converges for the weak topology (j{Aib, Co). 

Since we focus on scalar conservation laws, we can assume without loss of generality that 
the total mass of the system is scaled to 1 and thus we will work in some space of probability 
measures, namely the Wasserstein space of order g > 1, which is the space of probability 
measures with finite order q moment: 

Vq{M.^) = nonnegative Borel measure, /i(M^) = 1, J \x\'^ix{dx) < 00 

This space is endowed with the Wasserstein distance defined by (see e.g. |33[ IM] ) 

dwq{fJ',^)= inf { [ \y - x\'^j{dx,dy)\ (2.1) 
j&M [J J 

where r(/i, u) is the set of measures on x with marginals fi and u, i.e. 

r(/i,z/) = |7 G P,(M^ X M^); G CoiR''), J i{yMdy^,dy^) = j e(?/o)/^(c??/o) 

^{yih{dyo,dyi) = / ^{yi)u{dyi 



From the Kantorovich theorem, we know that in the definition of dwq the infimum is actually a 
minimum. A map that realizes the minimum in the definition (12. ip of dwq is called an optimal 
map, the set of which is denoted by To{n, u). 

A fundamental breakthrough in the use of the geometric approach to solve PDE is the work 
of F. Otto |28) . which is the basis of the so-called Otto Calculus jS^. Let X be a Riemannian 
manifold endowed with the Riemannian metric gx{-,-) (a positive quadratic form on the tangent 
space at X in X denoted T^.X). Let W : X — )■ M be differentiable. The gradient of W at x G X 
is defined by: for all v G T^X, let 7(t) be a regular curve on X such that 7(0) = x and 7'(0) = v, 
then 

^|t=oW(7)(t) = g,iVxW,v), V,.W G T,X. 
dt 

The gradient flow associated to W is the solution p : [0, +00) — )■ X of the differential equation : 

A fundamental result due to Ambrosio et al. p] states that gradient flows are equivalent to 
curves of maximal slope. Therefore, solving a PDE model of gradient type boils down to prove 
the existence of a curve of maximal slope. 

Let us be more precise. In the following, we will mainly focus on the case q = 2 and we will 
shortly denote dw instead of dw2- In the formalism of |2], we say that a curve fi is absolutely 
continuous, and we denote /i G AC^((0, +00), ^2(1^^)), if there exists m G L^(0, +00), such 
that div{lJ.{s), jj{t)) < J^m(r)dr, for < s < t < +00. Then we can define the metric derivative 

I r rfH/(/i(s),/i(t)) 

l/i \ {t) := limsup - 



s—^t 



s-t 
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The tangent space to a measure /i G V2(^ ) is defined by the closed vector subspace of L (fi) 



Tan^P2(K ) := {V0 : € C,°°(M^)} 

We recall the result of Theorem 8.3.1 of |2]: if /i G AC^{{0, +cx)), p2(ffi^)), then there exists 
a Borel vector field v{t) G L'^{fi{t)) such that 

dtH + V ■ (f yu) = 0, in the distributional sense on (0, +oo) x M^. (2.2) 

Conversely, if fi solves a continuity equation for some Borel velocity v G L^((0, +oo); L^(/i)) 
then fi is an absolutely continuous curve and |yu'|(t) < ||t;(t)||/^2(^). As a consequence, we have 

\fi'\(t) = min {||?;||z,2(^(i)) : v(t) satisfies ([22])} • 

Definition 2.1 (Gradient flows) We say that a map /i G y4Cf^^((0, +oo); 7^2(1^^)) a solu- 
tion of a gradient flow equation associated to the functional W if there exists a Borel vector 
field V such that v{t) G Tan^(^t)'P2{^^) for a.e. t > 0, ||i;(t) ||i2(^) G Lf^^{0, +00) , the continuity 
equation 

dtfi + V ■ (vfi) = 0, 

holds in the sense of distributions and v(t) G — 9W(p(t)) for a.e. t > 0, where 9W(p) is the 
suhdifferential ofW at the point p. 

We recall that the sub differential of a functional W is defined by 

w G dWifi) W{iy) - W(/i) > / w{x){e{x) - x) n{dx) + o(c/v^(/i, z/)). 



where 9 is an optimal transport map such that v = O^fi. Next, the slope |9VV| is defined by 

\R^Am ^ 1- (W(/i) - W(z/))+ 

\dW\ (/i) = hm sup — r , (2.3) 

where m+ = max{M,0}. We have the property ([2], Lemma 10.1.5) 

|9W|(/i) = min{||u;|U2(^) : w G 9>V(/i)}. (2.4) 

Moreover, there exists a unique w G VV(/i) which attains the minimum in (12. 4p . It is denoted 
by aoW(;u). 

Definition 2.2 (Curve of maximal slope) A curve fi G ACf^^{{0, +oo);V2(^)) is a curve of 
maximal slope for the functional Wift\-^ >V(/i(t)) is an absolutely continuous function and if 
for every < s < t < T , 

1 J'^ I/if (r) dr + y^ \dW\\f^ir))dr < W(/i(.)) - W(/i(t)). 

Finally Theorem 11.1.3 of |2] shows that curves of maximal slope with respect to {dWl 
are equivalent to gradient fiow solutions. Moreover, the tangent vector field v(t) is the unique 
element of minimal norm in the subdifferential of W (see fl2.4p ): 

v{t) = _aO>V(/x(t)) for a.e. t > 0. 
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2.2 Strategy of the proof in [14] 

The idea of the work by Carrillo et ah [H] is to extend the work of [2j to an interaction energy 
W defined through the interaction potential W in (ll.ip . whose derivative has a singularity in 
0. More precisely, attractive "pointy potentials" are considered, which we define now. 

Definition 2.3 (pointy potential) The interaction potential W is said to be an attractive 
pointy potential if it satisfies the following assumptions. 

(AO) W is continuous, W{x) = W{—x) and W{0) = 0. 

(Al) W is X-concave for some A > 0, i.e. W{x) — is concave. 

(A2) There exists a constant C > such that W{x) > — C(l + |xp) for all x G M^. 

(A3) W^gCHM^\{0}). 

Given a continuous potential IV : M — M, we define the interaction energy in one dimension by 

Mp) = \ I W{x-y) p{dx)p{dy). (2.5) 



2 _ 

The existence and uniqueness result of |14j can now be synthetized as follows. 



Theorem 2.4 (Theorem 2.12 and 2.13 of [T^.) Let W satisfies assumptions (AO)— (A3) and 

let a = id. Given p*™ G V2(^^), there exists a gradient flow solution of fll.ip . i.e. a curve 
p G ACL([0, oo);P2(ffi^)) satisfying 

^ + d,iv{t)pit)) = 0, zn V'i[0, oo) X R^), 
v{t,x) = -a°W(p)(t,x) = [ W'{x-y)p{t,dy), 

with p(0) = p*"*. Moreover, if pi and p2 are such gradient flow solutions, then there exists a 
constant A such that, for all t > 

dwiPi{t),P2{t)) < e^'dw{pi{0),P2{0)). 

Thus the gradient flow solution of fll.ip with initial data p*"* G 7^2 (IR^) '^s unique. Moreover, 
the following energy identity holds for all < to < ti < oo: 

\d^W * p\^p{t, dx)dt + W(p(ti)) = W(p(to))- (2.6) 

to 

Proof. We summarize here the main steps of the proof and refer the reader to for more 
details. The first step is to compute the element of minimal norm in the subdifferential of W. 
By extending Theorem 10.4.11 of [2], the authors prove that 

- 9°W(p) = d^W * p, (2.7) 
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where 



d°W*p{x)= f W'{x-y)p{dy). 



The second step is based on the so-called JKO scheme introduced in [23] (see also |2]). It 
consists in the following recursive construction for curves of maximal slope. Let r > be a small 
time step, we set = p*™ the initial data for (11. ip . Next, knowing p^, one proves (Proposition 
2.5 in [11]) that there exists p],j^^ such that 



p^+i G arg min j W(p) + ^(f^ipl, p) \ . {2.i 



Next, a piecewise constant interpolation p"^ is defined by 

p-(0)=p^™; p^{t)=pl if te{kT,{k + l)rl 

and Proposition 2.6 states the weak compacteness (in the narrow topology) of the sequence p"^ 
as r — 7- 0. Finally Theorem 2.8 ensures that the weak narrow limit p is a curve of maximal 
slope. 

The conclusion follows by applying Theorem 11.1.3 of [2], which allows therefore to get the 
existence of a gradient flow for the functional W. By definition, the gradient flow is a solution 
of a continuity equation whose velocity field is the element with minimal norm of the subdiffer- 
ential of W. In the first step of the proof this element has been identified to be d^W * p. Thus, 
it is a weak solution of the problem and moreover we have the energy estimate (12. 6p . □ 



2.3 The one-dimensional case 

We gather here several remarks specific to the one-dimensional framework. First we notice that 
assumptions (Al) and (A3) imply that x t— )■ W'{x) — Ax is a nonincreasing function on ]R\{0}. 
Therefore lim^.^o± W'{x) = W'{0^) exists and from (AO), we deduce that W'{0~) = -W'{0+). 
Moreover, for aA\ x > y in M\{0} we have W'{x) — Xx < W'{y) — Xy. Thus we have the one-sided 
Lipschitz estimate (OSL) for W 

Vx>yGM\{0}, W'{x)-W'{y) <X{x-y). (2.9) 

Letting ?/ -> 0^ we deduce that for all x > 0, W'{x) - Xx < W'{0+) and W'{x) - Xx < W'{0~). 
Thus we also have the one-sided estimate 

W'{x) < Xx, for all x > 0. (2.10) 

After an integration with (AO), we deduce that W{x) < Ax^/2. By concavity of the function 
W : X W{x) — Ax^/2, we have 

W'{x) h > W{x + h)- W[x) > -C{1 + x^ + h^), 

for any x 7^ 0. Taking h = max{|x|, 1}, we deduce that there exist nonnegative constants Co 
and Ci such that 

\W'{x)\ <Co + Ci\x\. (2.11) 
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The one- dimensional framework also allows to simplify several proofs in Theorem 12.41 In- 
deed any probability measure /i on the real line M can be described in term of its cumulative 
distribution function F{x) = oo,x)) which is a right-continuous and nondecreasing func- 
tion with F[—oo) = and F(+oo) = 1. Then we can define the generalized inverse of F by 
F~^[z) = inf{x G M./F{x) > z}, it is a right-continuous and nondecreasing function as well, 
defined on [0, 1]. If /i and u belongs to P2(K), with generalized inverses F~^ and G~^, we have 
this explicit expression of the Wasserstein distance (see [M] ) 

dw{fi,i^f= [ \F^\z) - G~\z)\^ dz. (2.12) 
Jo 

Moreover, we have 

{F''y{zUF-\z)) = l. 

In this framework, we can then rewrite the JKO scheme ( 12. Sp in the proof above. Let us 
denote by F^ the cumulative distribution of the measure pi and by Vj[ := (F^)~^ its generalized 
inverse. Then, in term of generalized inverses, (12.81) rewrites 

where ^ ^ 

W{V) =11 W{V{y) - V{z)) dydz. 
Jo Jo 

Such an approach using the generalized inverse has been used in [8] for the one dimensional 
Patlack-Keller-Segel equation. 



3 Duality solutions 

We turn now to the alternative notion of weak solution we wish to investigate. It is based 
on the so-called duality solutions which were introduced for linear advection equations with 
discontinuous coefficients in [9]. Compared with the gradient flow approach, this strategy allows 
a more straightforward PDE formulation, including from the numerical viewpoint. The main 
drawback is that presently duality solutions in any space dimension are only available for pure 
transport equations (see [H]). Since we have to deal here with conservative balance laws, we 
have to restrict ourselves to one space dimension. First we give a brief account of the theory 
developed in j^, summarizing the main theorems we shall use, next we define duality solutions 
for ([n]). 

3.1 Linear conservation equations 

We consider here conservation equations in the form 

dtp + d,{b{t,x)p) =0, {t,x) e{0,T) xR, (3.1) 

where 6 is a given bounded Borel function. Since no regularity is assumed for b, solutions to 
(13. ip eventually are measures in space. A convenient tool to handle this is the notion of duality 
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solutions, which are defined as weak solutions, the test functions being Lipschitz solutions to 
the backward linear transport equation 

dtp + bit, x)d^p = 0, p(T, .) = / G Lip(M). (3.2) 

In fact, a formal computation shows that ^ (^J^p{t,x)p{t,dx)^ = 0, which defines the duality 
solutions for suitable p's. 

It is now quite classical that to ensure existence for (13. 2p . it suffices for the velocity field to 
be compressive, in the following sense: 

Definition 3.1 We say that the function b satisfies the one-sided Lipschitz (OSL) condition if 

dxb{t, .) < /3{t) for P G L^{0,T) in the distributional sense. (3-3) 

However, to have uniqueness, we need to restrict ourselves to reversible solutions of (13. 2p : let 
£ denote the set of Lipschitz continuous solutions to (13. 2p . and define the set £ of exceptional 
solutions by 

6^ = |p G £ such that = o|. 
The possible loss of uniqueness corresponds to the case where S is not reduced to zero. 

Definition 3.2 We say that p & C is a reversible solution to liS. 2\] if p is locally constant on 
the set 

Ve = |(t,x) G [0,T] X M; 3 e ^, pe{t,x) ^ o}. 

We refer to [9] for complete statements of the characterization and properties of reversible 
solutions. Then, we can state the definition of duality solutions. 

Definition 3.3 We say that p G Sm '■= C{[0,T]; Mb{M.) - a{Mb,Co)) is a duality solution 

to ^3. 1\] if for any < t < T , and any reversible solution p to liS. 2\) with compact support in 

X, the function t t— / p{t,x)p{t,dx) is constant on [0,r]. 



We summarize now some properties of duality solutions that we shall need in the following. 

Theorem 3.4 (Bouchut, James J^) 

1. Given p° G A^fc(M), under the assumptions liS. 3\) . there exists a unique p G Sm, duality 
solution to lis. such that p{0, .) = p° . 

Moreover, if p° is nonnegative, then p(t, ■) is nonnegative for a.e. t> 0. And we have the 
mass conservation 

|p(t,-)|(K) = |p1(ffi), for a.e. t e]0,T[. 

2. Backward flow and push-forward: the duality solution satisfies 

Vt G [O,T],V0 G Co(M), [ ^{x)pit,dx) = [ (j){X{t,0,x))p%dx), (3.4) 

where the backward fiow X is defined as the unique reversible solution to 
dtX + b{t, x)dxX = zn]0,s[xM, X{s,s,x)=x. 
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3. For any duality solution p, we define the generalized flux corresponding to p by h^p = 
—dtu, where u = pdx. 

There exists a bounded Borel function b, called universal representative of b, such that 
b = a almost everywhere, and for any duality solution p, 

dtp + dx{bp) = in the distributional sense. 

4. Let ibn) be a bounded sequence in L°°(]0, T[xM), such that bn ^ b in L°°(]0, T[x]R) — wi^. 
Assume dj)n < ttn(^); where is bounded in L^(]0,T[), d^b < a E L^(]0,T[). Consider 
a sequence (p,„) G Sj^ of duality solutions to 

dtPn + dx(bnPn) = ^ m ]0,T[xM, 

such that pn{0, .) is bounded in and p.„(0, .) p° G A^b(M). 

Then p^ ^ p in Sm, where p G Sm is the duality solution to 

dtp + dx{bp) = Q m ]0,r[xM, p(0,.) = P°. 

Moreover, bnPn bp weakly m A^fc(]0, T[xM). 

The set of duality solutions is clearly a vector space, but it has to be noted that a duality 
solution is not a priori defined as a solution in the sense of distributions. However, assuming 
that the coefficient b is piecewise continuous, we have the following equivalence result: 

Theorem 3.5 Let us assume that in addition to the OSL condition (13. 3p . b is piecewise con- 
tinuous on ]0,T[xM where the set of discontinuity is locally finite. Then there exists a function 
b which coincides with b on the set of continuity of b. 

With this b, p E Sm is a duality solution to ( 13. ip if and only if dtp + dx{bp) = in P'(M). 
Then the generalized flux b^p = bp. In particular, b is a universal representative ofb. 

This result comes from the uniqueness of solutions to the Cauchy problem for both kinds of 
solutions (see Theorem 4.3.7 of [H]). 

3.2 Duality solutions for aggregation 

Equipped with this notion of solutions, we can now define duality solutions for the aggregation 
equation. The idea was introduced in [TU] in the context of pressureless gases. It was next 
applied to chemotaxis in [21], and we shall actually follow these steps. 

Deflnition 3.6 We say that p G C{[0,T]; Aib{^)) is a duality solution to (II. ip if there exists 
hp G L°°((0,T) X M) and a G L;^„^(0,T) satisfying d^cip < a in V'{{0,T) x R), such that for all 
<ti<t2<T, 

dtp + dx{dpp) = in the sense of duality on (^1,^2), (3-5) 
and CLp = a{W' * p) a.e. 
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This allows at first to give a meaning to the notion of distributional solution, but it turns out 
that uniqueness is a crucial issue. For that, a key point is a precise definition of the product 
OpP, as we shall see in more details in Section 13.31 below. 

We now state the main theorems about duality solutions for the aggregation equation f ll.ip . 
Existence of such solutions in a measure space has been obtained in [21] in the particular case 
W{x) = ^e~'^' and a similar result is presented in |19| when W{x) = —\x\/2 which appears in 
many applications in physics or biology. We extend here these results for a general potential 
satisfying assumptions (AO)— (A3). However to do so, we have, as in |[T4], to restrict ourself to 
the linear case, that is a = id. 

Theorem 3.7 (Duality solutions, the linear case) Let W satisfies assumptions (AO)— (A3) 
and let a = id. Assume that p*™ G Pi(]R). Then for any T > 0, there exists a unique duality 
solution p>0 of equation (II. ip in the sense of Definition \3. 61 with p(0) = p*™, p{t) G Pi(M) for 
any t G (0, T) and which satisfies (13. 5p in the distributional sense with 

ap{t,x) := d'^W * p{t,x) = W'{x -y)p{t,dy). (3.6) 

Jx^y 

Then we have p = X^p*"* where X is the backward flow corresponding to cip. 

We turn now to the case a ^ id. In order to be in the attractive case, we assume that the 
function a satisfies the following 

Assumption 3.8 a is non- decreasing, with 

aGC^(M), 0<a'<a. (3.7) 

In this context, existence and uniqueness of duality solutions have been proved for the case of 
an interaction potential W = |e~l^'l in [21]. We extend here the techniques developed in this 
latter work to more general potentials W. However we are not able to prove such results in 
the whole generality of assumptions (AO)— (A3) and need more regularity on the interaction 
potential. 

Assumption 3.9 Let W ^ C^{W \ {0}), we assume that in the distributional sense 

W" = -5o + w, we Cfe°(M), (3.8) 

where 6o is the Dirac measure in 0. 

This allows a definition of the flux in (11. ip which generalizes the one in [21]. Indeed we can 
formally take the convolution of (13. 8p by p, then multiply by a{W' * p). Denoting by A the 
antiderivative of a such that A(0) = and using the chain rule we obtain formally 

- d^AiW' *p) = -a{W' * p)W" *p = a{W' * p){p + w*p). (3.9) 

Thus a natural formulation for the flux J is given by 

J := -d^A{W' * p) + a{W' * p)w * p. (3.10) 
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Theorem 3.10 (duality solutions, the nonlinear case) Let us assume that is given in 
'Pi(R). Under Assumption \3.9[ on the potential W and (13. 7p for the nonlinear function a, for all 
T > there exists a unique duality solution p of (II. ip in the sense of Definition \3.6[ p(t) G Pi(M) 
for t G (0,T) and which satisfies in the distributional sense 

dtp + d,J = 0, (3.11) 

where J is defined by (I3.10p . 

Theorems 13.71 and 13.101 are proved respectively in Sections H] and [5l but before diving into the 
detailed proofs, we comment the main steps, which are common to both cases. 

• existence of duality solutions is obtained by approximation. First we obtain the dynamics 
of aggregates (that is combinations of Dirac masses), then we proceed using the stability 
of duality solutions 

• uniqueness is obtained by a contraction argument in Pi(R). No uniqueness is expected 
in a general space of measures. The argument is repeated in V2{^) so that gradient flow 
and duality solutions can be compared. 

3.3 Velocities and fluxes 

When concentrations occur in conservation equations, leading to measure-valued solutions, a 
key point to obtain existence and uniqueness in the sense of distributions is the definition of 
the flux and the corresponding velocity. This was already pointed out in [10], where duality 
solutions are defined for pressureless gases, and partially managed through conditions on the 
initial data. A more satisfactory solution came out in since uniqueness was completely 
handled by a careful definition of the flux of the equation, or in other terms, the product dpp. 
An analogous situation arising in plasma physics is considered in [19], around duality solutions 
as well. In a similar context, other definitions of the product can be found, see |27| in the 
one- dimensional setting, and [32] for a generalization in two space dimensions, where defect 
measures are used. 

We explain in more details this point in our context, in order to give a meaning to both 
duality and gradient flow solutions in the sense of distributions. As a rule, the product of 
a{W' * p{t)) by p{t) is not well-defined when p{t) G A4b(M). First we compute W * p. We write 
p = dxU, so that u G BV{W). For such a function, we denote by Su the set of a: G M where u 
does not admit an approximate limit, \Su\ = 0, and by J„ C Su the set of approximate jump 
points (see Remark 3.98 of [T]). We use the decomposition p = d^u + p'^ + fP , where « C 
is the regular part of the derivative, fP = 'Yliy^j Cy^y ^^e jump part, and p'^ the so-called Cantor 
part. The so-called diffuse part of the derivative is p'^ = d^u + p^. For x ^ J„, we easily obtain 

W' * p{x) = W'* p\x) + CyW\x - y), 

y&Ju 

while if X G J^, the function is not defined. Indeed, letting — )■ x, first with z < x, then with 
z > X, we obtain 

W' *p{x^) = W' *dlu{x)+ CyW\x-y) + CxW'{Q^). (3.12) 
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Removing the indetermination amounts to define a velocity for a single Dirac mass located in 
X or equivalently for the center of mass of the density. Obviously, formula (13. 6p sets this value 
to 0, hence a single Dirac mass is stationary, and the product by the measure p is meaningful. 
Therefore in the linear case we can consider the fiux 

J(t, x) := ap{t, x) pit, x). 

Recall that this value is obtained by computing the element of minimal norm in the subgradient 
of the energy W corresponding to W . 

On the other hand, in the nonlinear case, with W" = —6q + w, the natural quantity to be 
defined is the fiux J, by formula (I3.10p . To define the corresponding velocity, and give rigorous 
meaning to (13. 9p . we use the Vol'pert calculus for BV functions |35] (see also [1]): for a BV 
function u, the fonction ay defining the chain rule dxA{u) = ayd^u is constructed by 

avix) = I a{tui{x) + [1 — t)u2{x)) dt, (3.13) 

where 

{(m, u) if X G M \ Su, 

{u^,u~) if X G Ju-, (3-14) 

arbitrary elsewhere. 
Now we apply that to u = W * p, and obtain, using the antiderivative A of a, 

' a(iy' *p(x)) if X G M\5«, 



ay(x) 



A{W' . p{x^)) - A{W' . p{x~)) .3^5. 
arbitrary elsewhere. 



The connection with the linear case follows since then Aiv) = hence 
A{W' * p(x+)) - A{W' * p(x-)) W * p(x+) + W'* p{x- 



W * p{x+) - W * p{x~) 2 

Therefore the undetermined term in (I3.12p is replaced by {W'{0^) + W'{0~))/2, which vanishes 
since W is even, and we recover (13. 6p . 



4 The linear case 

By linear case we mean the case where a = id in (II. ip . Together with assumptions of Definition 
12. 3[ this is exactly the context of [H]. At first we prove Theorem 13. 7[ obtaining existence of 
duality solutions in Subsection 14. H and uniqueness in Subsection 14.21 Next, in Subsection 14.31 
we establish that they are equivalent to gradient fiow solutions, thus answering the questions 
raised by Remark 2.16 of |14j . More precisely, we prove the following theorem. 

Theorem 4.1 Let a = id. Let us assume that W satisfies assumptions (AO)— (A3) and that 

pini g p2(M). 
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(i) Let p be the duality solution as in Theorem \3. 71 Then for all t > 0, p(t) G V2(^), 
p G AC f^^{{0, +00); V2{^)) and p is the gradient flow solution as in Theorem^KJ. 



(a) If p is the gradient flow solution of Theorem 2.4 , then it is a duality solution as in Theorem 
4.1 Existence of duality solutions 

The first step is to verify that the velocity ap defined by (13.61) satisfies the OSL condition (13. 3p . 

Lemma 4.2 Let p{t) G A^f,(M) be nonnegative for all t > 0. Then under assumptions (AO) — 
(A3) the function {t,x) 1— )■ dp{t,x) defined in (13. 6p satisfies the one-sided Lipschitz estimate 

ap{t,x) - dp{t,y) < X{x - y)\p\(R), for all x>y,t>0 

Proof. By definition (13. 6p . we have 



dp{x) - dp{y) = / {W\x -z)- W\y - z))p{dz) + W\x - y) / p{dz), 

J z^x,zi^y J z={x}\j{y} 

where we use the oddness of W (AO) in the last term. Let us assume that x > y, from (12. 9p . we 
deduce that W'{x — z) — W'{y — z) < \{x — y) and with (I2.10p . we deduce W\x — y)< \{x — y). 
Thus, using the nonnegativity of p, we deduce the one-sided Lipschitz (OSL) estimate for dp. □ 

Proof of the existence result in Theorem 13.71 This proof is split in several steps. 
• Aggregates. 

Let us assume in a first step that p*"* = X]r=i '^i^x'^ where x\ < x\ < ■ ■ ■ < x^^ and the mj-s 
are nonnegative. We look forward a solution pn{t,x) = Yyi=i'^i^x,{t) in the distributional sense 
of the equation 

dtp + d^idpp) = 0, (4.16) 

where dp is defined in (13. 6p . Let n„ := Pn = Yl^=i''^i^i^ ~^ ^iit)) where H is the Heaviside 
function. Then we have that 

n 

-dtUn = dp^pn = ^mi^mjW'{xi - Xj)S^^. 



In fact, 



^^^mjW'{xi — Xj) ii X = Xi, i = 1 



n 



mjW'{x — Xj) otherwise. 

i=i 

Then the sequence {xi)i=i^,,,^n satisfies the ODE system 



(4.17) 



x'iit) = ^mjW'{xi-Xj), Xi(0) = x°, i = l,...,n£, (4.18) 
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where rii < n is the number of distinct particles, i.e. rii = ^{i G {1, . . . , N},Xi ^ Xj,'^j}. 
Notice that fl3.12p rewrites as 



W' * pn{xf) = ^TTijW'ixi - x^) + miiy'(0±). 
We define the dynamics of aggregates as follows. 

• the Xj-s are solutions of system ( I4.18P (where the right-hand side is zero if ni = 1), when 
they are all distinct; 

• at collisions, we define a sticky dynamics: if Xi = Xj at time Ti when for instance i < j, 
then the two aggregates collapse in a single one and we redefine system fl4.18p by changing 
the number to — 1, replacing the mass rrii by rrii + rrij and deleting the point Xj. We 
denote by := Tq < Ti < . . . T^. < oo the times of collapse, where k < n. 

This choice for the dynamics is clearly mass-conservative. Moreover, with this choice at 
the coUisional times, we still have that Xj < Xj when i < j. Thus the function defining the 
right-hand side satisfies the one-sided Lipschitz condition and there exists a solution to this 
dynamical system. 

In fact, the function defining the right-hand side of fl4.18p is continuous and locally bounded 
as long as all Xi are distincts. Therefore there exists a local solution of the ODE f l4.18p 
thanks to the Cauchy-Peano Theorem. A difficulty appears when particles collide. Thanks to 
the one-sided estimate on W (12. 9p . we have that the right hand side defining fl4.18p is one- 
sided Lipschitz. Hence, there exists a unique Filippov fiow which is global in time thanks to 
estimate (12. lip . Moreover this unique solution coincides with the solution obtained by the 
Cauchy-Peano theorem on each time interval [Ti,Ti^i). Finally we have constructed a solution 
{xi)i=i^,,,^n of the system of ODE (I4.18P which is Lipschitz on [0,T] and on each interval 
{Tt, Ti+i). 

Thus we can define p„ := Y2^=i''^i^xi- It is then straightforward by contruction that p is a 
solution in the distributional sense of (I4.16p - (l3.6p . 

• Duality solutions. 

Using Lemma [4.21 and expression (I4.17p . we deduce that cip^ satisfies the OSL condition and 
is piecewise continuous where the set of all discontinuities is {xj}j=i ... From Theorem 13. 5[ we 
deduce that p„ is a duality solution for all n G N*. 

• Finite first order moment. 

By construction, we clearly have that p > 0. Let us denote by ji{t) := Yl^=i''^i\^ii^)\ 
first order moment. We have ji(0) < +oo. Using (I4.18p . we compute 

n n 

j'i{t) < ^mi^mj\W'{xi -Xj)\ < C^mi^mj(l + \xi - Xj\), 

i=l j^i i=l jjki 

where we use (12. lip for the last inequality. Using 'Y^^mi = |p*"*|(M) = 1, we deduce that there 
exists a nonnegative constant C such that 

m<c{i+m). 
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Integrating this latter inequality implies that ji is bounded on [0, T] by a constant only depend- 
ing on T and on ji(0). It remains to pass to the limit n — +00 in the regularization. 
• Passing to the limit n — )■ +00. 

Using the fact that Pri(^) G ^i(I^) ior all t > and from f l2.1ip . we deduce that cip^ is 
bounded in L°°((0, T) x M) uniformly with respect to n. Thus, from the point 4 of Theorem 13 .4^ 
we can extract a subsequence a„ converging in L°°((0, T) x M) towards a and the corresponding 
duahty solutions sequence {pn)n converge in Sm towards p which is a duality solution of 

dtp + d^{ap) = 0. 

Moreover, since p„ ^ p in Sm , we have that dp„ — ?■ Op a.e. with the definition in (13.6p . Thus 
a = ttp a.e. and the flux 

J„(t, x) J{t, x) := hpp in yWbQO, T[xM) - w;. 
4.2 Uniqueness 

Let p be a nonnegative duality solution which satisfies in the distributional sense 

dtp + d^{(Xpp) = Q. (4.19) 

As above, we denote by F the cumulative distribution function of p and by F~^ its generalized 
inverse. We have then by integration of fl4.19p 

dtF + a,AF = 0. 

Then, we deduce from standard computations that 

dtF~\t, z)p{t, F-\t, z)) = apit, F-\z))p{t, F-\t, z)). 

Then the generalized inverse is a solution of the system 

dtF-\t,z) = dp{t,F-\z)). (4.20) 

Moreover thanks to a change of variables in f l3.6p . 

dp{t, F-\z)) = I W\F-\z) - F~\y)) dy. 

Proposition 4.3 Let pi(t, ■) and P2{t,-) belong to 7^1(1^). If pi and p2 satisfy equation fl4.19p 
in the weak sense with initial data p^^"* and p^^\ Then we have, for all t > 

dwi{pi{t, ■),p2{t, ■)) < e^^dwi{pup2)- 

Proof. Let us assume that there exists two nonnegative measures pi and p2 duality solutions 
which satisfy moreover fl4.19p . Let us denote respectively F^^ and ^2"^ their generalized inverses. 
From fl4.20p . we have 

dt{F,-' - F,-') = dp{t,F,-\z))-dp{t,F,-\z)). 
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Multiplying the latter equation by sign[F^ {z) — F2 (z)) and integrating, we get 
d,{Fr\z)-F,-\z))sign{Fr\z)-F,~\z)) = {d,{Fr\z))-a,{F,^^^^^^ 

(4.21) 

The one-sided Lipschitz estimate on Sp in Lemma [4.21 implies that 

{a,{Fr\z)) - a,{F,-\z)))sign{Fr\z) - F^\z)) < X\F,~\z) - F,-\z)\. 
Hence, after an integration, 

1 1' \F,~' - F,-'\{z) dz<\f^ iFfi - F^\z) dz. 
Applying a Gronwall argument, we deduce that 

\\{F{' - F^'){t)\Wi^,i) < e''m~' - F2"^)(0)|k.(o,i) 

We conclude the proof by noticing that || {F{'^ — -F2~^)(t) ||li(o,i) = dwi{pi, P2) and using a Gron- 
wall argument. □ 

Proof of Theorem 13.71 The existence has been obtained in Section 14.11 Then if we have 
two duality solutions pi and p2 as in Theorem 13.71 Proposition 14 . 3 1 implies that their generalized 
inverse are equals. Therefore pi = p2- Finally, the second point of Theorem 13.41 allows to define 
the duality solution as a push-forward measure from the backward flow. |--| 



4.3 Proof of Theorem 14.11 

Let us first prove that the second order moment is bounded provided p*"* G V2{9.). We follow the 
idea of the proof of finite first order moment in subsection 14.11 : we consider an approximation of 
pim Yl^=i '^i^x° and we consider then the constructed duality solution p(t, x) = Y17=i "^i^Xiit) 
where the dynamics of the nodes {xj}j=i^...^„ is given in fl4.18p . Let us denote by j2{t) '■= 
the second order moment. We have j2(0) < +00. Using f l4.18p . we compute 

n n 

i2{t) = 2mjXj mjW'{xi — Xj) < C "2j(l + \xi — 

i=l j^i 1=1 jj^i 

where we use fl2.1ip for the last inequality. Since ^^mj = Ip^^KM) and from the Cauchy- 
Schwarz inequality, we deduce that there exists a nonnegative constant C such that 

n 

m < C5^m,(|x,| + |x,n < C(l + 
1=1 

Integrating this latter inequality implies that j2 is bounded on [0,T] by a constant only de- 
pending on T and on J2(0). Then we can pass to the limit n — )■ +00 to obtain a bound on 
Jjj |a;pp(t, dx) for any t > 0. Moreover, using this L°°((0, T); ^2(1^)) bound, we deduce that the 
velocity field dp defined in (13.60 is bounded in L^((0,T); L'^{p)). 
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Then, if p is a duality solution as in Theorem \3.7\ it satisfies in the distributional sense the 
scalar conservation laws 

dtp + d^{app) = 0. 

Using Theorem 8.3.1 of |2] and the L^((0,T); L'^{p)) bound on the velocity hp, we deduce that 
p G ACf^^{{0, +oo);V2{^)) ■ Thus p is a gradient fiow solution (see jT^ or Sections 8.3 and 8.4 
of |2]). It concludes the proof of point (i) of Theorem 14.11 

Conversely, if p is a gradient fiow solution of Theorem 12. 4[ then it is a weak measure solution 
of f l3.5p -f l3l6|) . Such solution is unique from Theorem 13.71 then p is a duality solution. □ 



5 On the case a ^ id 

In many applications in physics and biology (see e.g. [TII| 1^). we have to consider that a ^ id. 
This situation is not so favourable as the previous one, and one has to impose restrictions on 
the potential W. First we recall that attractivity implies that a is non-decreasing, see (13. 7p . 
Next, we need to assume that W has the following structure, see (13. 8p . 

W" = -6o + w, we C° 

With these assumptions, we are able to prove existence and uniqueness of duality solutions. 
Theorem I3.10[ this is the aim of subsection 15.11 Next, in subsection 15. 2[ we turn to gradient 
fiow, which are definitely not well suited for that case, since we have to restrict ourselves to 
w = in the previous assumption on W. 

5.1 Duality solutions 

Here we prove Theorem 13.10^ following the same strategy as in the linear case: first we prove 
the OSL condition, next establish the dynamics of aggregates, which leads to existence by 
approximation. Finally, uniqueness follows from a contraction principle in the space Vi. In 
addition, we prove that duality solutions are absolutely continuous in time. 

5.1.1 OSL condition 

The first step consists in checking the OSL property for a. 

Lemma 5.1 Assume < p G Aib(^) o-nd that (13. 7p holds, If Assumption [XPI is satisfied, then 
the function x a{W' * p) satisfies the OSL condition (13. 3p . 



Proof. Using (13. 8p . we deduce that 

dxxW * p = -p + w * p. 

Therefore, 

dx{a{dxW * p)) = a'{dxW * p){—p + w * p) < a'{d.jW * p)w * p, 
where we use the nonnegativity of p in the last inequality. Then from (13. 7p we get 

dx{a{d,M * p)) < a\\p\\L^\\w\\L^. 
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It implies the OSL condition on the velocity. □ 

5.1.2 Existence 

Approximation with aggregates. 

Following the idea in subsection 14.11 we first approximate the initial data by a finite sum 
of Dirac masses: = XliLi ^i^x" where < < ■ ■ ■ < 2;° and the mj-s are nonnegative. We 
look for a sequence (pn)n solving in the distributional sense dtPn + dxJn = where the fiux J„ 
is given by (I3.10p . A function pn(t,x) = ^"=1 ?Tii5a;j(t) is such a solution provided the function 
Un defined by 

Un{t,x):= pndx = ^miH{x - Xi{t)), (5.1) 
^ 1=1 

where H denotes the Heaviside function, is a distributional solution to 

dtUn - dxA{dxW * pn) + a{dxW * pn)w * p„ = 0. (5.2) 

We have 

n 

W *pn{t,x) = ^m,W'{x-x,{t)). (5.3) 

i=l 

And from (13. Sp . we deduce that 

/X 
w{y)dy. (5.4) 

Straightforward computations show that in the distributional sense 

n 

dxA{dxW * Pn) = a{dxW * pn)w * Pn 

i=l 

where [/jj.- = f{xf) — f{x^) is the jump of the function / at Xj. Injecting (15. ip and (15. 5p in 
(IE2D, we find 

n n 

-^miXi(t)4,(i) = ^[A{dxW * pn)]x,5x^. 
1=1 1=1 

Thus it is a solution if we have 

mix[{t) = -[A{dxW * pn)]x,{t), for 2 = 1, . . . ,n. (5.6) 
This system of ODEs is complemented by the initial data Xj(0) = Moreover, we have from 

i n 

W' * Pn{xf) = — ^^"^j + ^^mju;(xj — Xj). (5.7) 
i=i j=i 

And 

W'*Pnix-)=m, + W'*pn{xl). (5.8) 
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Thus 

Then we define the dynamics of aggregates as in subsection 14.11 : 

• When the Xi are all distinct, they are solutions of system fl5.6p or equivalently fl5.9p (with 
zero right hand side if ne = 1). 

• At collisions, we use a sticky dynamics as defined as above. 

We recall that this choice of the dynamics allows the conservation of the mass. As above, we 
have existence of the sequence satisfying ( 15. 9p on [0,T] with initial condition Then 
we set pn(t,x) = J27=i''^i^Mt)(-^)- construction, p„ is a solution in the sense of distribution 
of f l3lT]) -f l3:TU|) for given initial data 
Existence of duality solutions 

We recall the following result due to Vol'pert [35] (see also [T]): if u belongs to i?V"(]R) and 
/ e C^{R) with /(O) = 0, then v = fou belongs to BV{R) and 

^lu with 7„ = f'{u) a.e. such that (/ o u)' = jy. 

Together with the fact that A is an antiderivative of a such that v4(0) = 0, this result implies 
that there exists a function a„ such that a„ = a{W' * pn) a.e. and 

dx{A{W' * Pn)) = ttndxxW * Pn = -(InPn + ttnW * p„, 

where the last identity comes from (13.81) . Therefore, we have 

Jn-.= -dx{A{W' * Pn)) + a{W' * pn)w * pn = anPn, and an = a{W' * Pn) Si.e. 
Thus Pn is a solution in the distributional sense of 

dtPn + dxittnPn) = 0. 

Moreover, we deduce from (15. 3p that a(W' * pn) is piecewise continuous with the discontinuity 
lines defined by We can apply Theorem 13.51 which gives that p„ is a 

duality solution and that a„ is a universal representative of a{W' * pn)- Then the fiux is given 
by a{W' * pn)^Pn = Jn- 
General case. 

Let us yet consider the case of any initial data p*"* G A^fe(M). We approximate p*™ by 
Pn^ — ^l=i''^i^x° with p™* ^ p*"* in A^fe(M). By the same token as above, we can construct 
a sequence of solutions (p„)„, with p„(t = 0) = p^"* = Yl^=i'^i^x°i which solves in the sense of 
distributions 

dtpn + dxJn = 0, Jn = -d^A^d^W * pn) + a^d^W * pn)w * Pn, 

and which satisfies 

a-npn = Jn, O-n = a{W' * p„) a.C. 

Moreover, since W * pn is bounded in L°° uniformly with respect to n by construction, we can 
extract a subsequence of {a{W' * pn))n that converges in L°° — weak* towards b. Since from 
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Lemma [5. H a{W' * pn) satisfies the OSL condition, we deduce from Tlieorem l3.4l 4) tliat, up to 
an extraction, p„ ^ p in Sm and a„p„ ^ ap weakly in A^;,(]0, T[x]R), p being a duality solution 
of the scalar conservation law with coefficient h. Then Jn ^ J := —dxA{W' * p) + a(W' * p)w * p 
in X''(M) and that a{W' * Pn) — ?■ a{W' * p) a.e. By uniqueness of the weak limit, we have 
h = a{W' * p). Moreover J = ap a.e. and p satisfies then fl3.1ip . 
Finite first order moment. 

As in subsection 14. we define ji{t) = XliLi compute 



1=1 



\Xi\ [W'*Pr, 



where we use (15. 9p . Since a is nondecreasing, A is a convex function. Moreover, we deduce from 
(15. 4p that < C(l + Thus, using the fact that a' is bounded (13. 7p . we have 

n 

<C7^m,(l + ^m,|x,-x,|) <C(l+ji(t)), 

i=l j 

where C stands for a generic nonnegative constant. Applying a Gronwall inequality and passing 
to the limit n +oo allows to conclude the proof. |--| 



Remark 5.2 Let us consider the case studied in the previous Section : a = id and W is even. 
Since W is odd, then (15. 4p rewrites 

r 1 

W' [x] = —H{x) + wo{x), where wq{x) = / w{y)dy+-. 

Jo 2 

When a = id, we have A{x) = Then system (15. 6p rewrites 

m,x[{t) = ~{W' * P„(X+) - W * Pnix-))iW' * Pnixf) + W * p„(x-)). 

Then, from { \5.7h and (15. Sp . we have 

= - ^^J - Y "^iW^o(a;i - Xj) = - ^mj + rrijWoixi - Xj). 

j = l j = l j = l jjLi 

From the expression of W above, we deduce that x[{t) = ^j^ifTijW{xi — xj) and we recover 
the dynamical system (I4.18P of the previous Section. 

Remark 5.3 The dynamical system (15. 6p defines actually the macroscopic velocity. Indeed, if 
we formally take the limit n — +oo of the right hand side of (15. 9p . this latter term converges 
towards the velocity ay defined by the chain rule (I3.15P 

5.1.3 Uniqueness. 

The proof of the uniqueness follows straightforwardly the one in subsection 14.21 since the macro- 
scopic velocity still satisfies the OSL condition as proved in Lemma [5. II 
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5.1.4 Absolute continuity. 



Following subsection 14. 3[ we have the following result: 

Proposition 5.4 Under the assumptions of Theorem \3.1C\ if moreover p™ G V2(^), then the 
duality solution of Theorem \3.1Ch satisfies p G ACf^^{{0, +oo); 



Proof. The proof of the finite second order moment follows straightforwardly the one for the 
first order moment in subsection 15.1.21 Then from (13. 7p and (12. lip , we have that 



\a{W'*p)\<c{l + Jjx-y\pidyi 
Applying a Cauchy-Schwarz inequality, we deduce that 

\a{W' * p)]"^ p{dx) < C (^1 + I \x - y]"^ p{dx) p{dy)y. 



which is bounded since p{t) G V2{^)- Therefore the velocity field a is bounded in L^{{0, T); L'^{p)). 
Moreover p is a solution in the distributional sense of dtp + dr^{ap) = 0. We conclude then from 
Theorem 8.3.1 of [2] that p G ACl^{{0, +00); P2(ffi))- □ 



5.2 Gradient flows 

When a 7^ id, equation (II. ip does not have a structure of a gradient. Moreover, we are not able 
to determine a conserved energy corresponding to the system. Nevertheless, in the particular 
case W{x) = —^\x\, we are able to adapt the technique of [H] to recover the existence of 
gradient flow solutions. We introduce for p G 7^2(1^) the functional 

W{p) = — f xa{u{x)) p{dx) , where u{x) = W' * p{x) = f p{dy) — -. (5.10) 
Jm. J-00 2 

Recall that we denote by A the antiderivative of a which vanishes at 0. We first verify that when 
a = id this energy is equal to the one introduced in (12. 5p . In fact, we have for W{x) = — 



-7/ \x-y\p{dx)p{dy) = -\ f [ {x-y)p{dy)p{dx) + \ f [ {x-y)p{dy)p{dx). (5.11) 
Repeated use of Fubini's Theorem in the last term of the right hand side leads to 
-7 / \x- y\p{dx)p{dy) =-]-!! {x - y)p{dy)p{dx) 

X 



\ I ^ I Pidy)p{dx) + ^ / / p{dx)yp{dy). 




Since u is given by (15.10p . we have J^°° p{dy) = | — u{x). Finally, 

— - I \x — y\p{dy)p{dx) = — f xu{x)p{dx 

which concludes the proof. 

We are able to prove in this case the following Theorem. 



22 



Theorem 5.5 Let W{x) = —^\x\ and a satisfy assumption {\3.7\i . Given G V2(^), there 
exists a unique gradient flow solution p G AC f^^{[0, +00), V2{^)) in the sense of Definition \2.1[ 
Therefore p satisfies in the distributional sense 

dtp + d^{a{u)p) = 0, with p(0) = p™, 

where u{x) = W * p{x) = p((— oo,x)) — |. Moreover, this solution is unique and we have the 
energy estimate : for all < to < ti < 00, 

ti 

2 



a{u{t,x))\' p{t,dx) dt + W{p{ti)) = W(p(to))- 



to 



Proof. This proof closely follows the ideas of the Section 2 of [13]. Therefore we will only 
gives the main steps and do not detail all the computations. As recalled in the beginning of 
this paper, the proof is divided into several steps. 

• —a{u) is the unique element of minimal L^(p)-norm in the subdifferential dW 
of W. 

We first show that —a{u) G 9W(p), where u is defined in flS.lOp . That means 



>V(%p)-W(p) > - / a{u){x){e{x) ~ x) p{dx) + o{dw2{0#P, p)), (5.12) 



where 9 is an optimal transportation map. From Brenier's Theorem [T^, we can take for 9 a 
nondecreasing function. By definition of the pushforward measure, we have 



r di9#p) = pi9-\-oo,x)), 
J —00 



so that 



yV{9#p) = - [ xa{p{9-\-oo,x)))9#p{dx) = - [ 9{x)a{p{{y G R/9{y) < 9{x)}))p{dx). 
We deduce from the monotony of the fucntion t that 

>V(%p) - >V(p) = - / a{u{x)){9{x)-x)p{dx). 



To prove that —a{u) is an element of minimal norm in 9VV(p), we first notice that for ^ G 
and for e > small enough such that {id + eC,) is increasing, we have 



00 / 



yV{{id + e^)#p) = - {id + e^){x)a{u{x))p{dx). 
Jr 

Thus, 

W{{id + eO#p) -^^{p) f , , ^ , , N 

lim — = - / a{u{x))i{x)dp{x). 

^-^0 ^ Jr 

Then, from the definition of the slope (12.31) . we have 

e\o dw2{{id + eQ#p, p) 
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We deduce from above that 

a{u{xmx) dp{x) < |9W|(p)liminf ^^^^:aM±£ik^ < |9W|(p)||e|U^(p), 

where we use the inequahty dw2{S#p,T^p) < — T||i2(p) for the last inequahty. By the same 
token with —C, instead of ^, we deduce 

aiuixmx) dpix)\ < IdWlipMhHp). 

Since ^ is arbitrary, we get ||a(M) ||j^2(p) < |9VV|(p). We conclude by using (12. 4p . 
• Well-posedness and convergence of the JKO scheme. 

The JKO scheme is defined in (12. 8p where r > is a given small time step. The convergence 
of this scheme can be deduced by an adaptation of the argument of Subsection 2.2 of |14] . 
However, in this one dimensional case, we can simplify a lot the computations by using the 
generalized inverse, denoted v, of the function u, in the spirit of [8]. In fact, the function u is, 
up to a constant, the cumulative distribution of p. Then, we have 

W{p) = W{v) := - / v{t, z)a{z) dz, 

J -1/2 

where we recall 

v{t, z) = u~^{t, z) := inf {x G R/u{t, x) > z}. 

We define moreover w*™ = (m*™)~-^, where m*™(x) = p*™((— oo, x)). Thus, using moreover fl2.12p . 
the minimization problem defined in (12.80 is equivalent to : let t;^ G L^(— 1/2, 1/2), we define 

vl, 1 = arg min ly^{v) + -^\\v - vlf] . (5.13) 

We recall that a is a nondecreasing function, then a(— 1/2) < a{z) < a(l/2) when z G 
(—1/2, 1/2). Then the functional defined inside the brackets is clearly lower semi-continuous, 
convex and coercive on L^(— 1/2, 1/2). We deduce that the above minimization problem (I5.13P 



admits an unique solution. Moreover, computing the Frechet derivative of the functional defining 
(I5.13P in the L^-metric, the Euler-Lagrange equation associated to this minimization problem 
implies 

1 1 

vl^^{z) =vl{z) -Ta{z), ^^^^2' 2^' 

This is an implicit Euler discretization of the equation dtv{t, z) + a{z) = 0. It is well 
known that the piecewise constant interpolation v'^ defined by v'^{0) = f*™ and v'^{t) = vj, if 
t G (/cr, {k + l)r], converges in L^(— 1/2, 1/2) as r — t- towards v{t) for all t G [0, +oo), where 
v{t, z) = w*™ — ta{z). Moreover, we have the energy estimate 

W{v) = W{ir')- / \a{z)\^dz. (5.14) 

Jo J~l/2 

We define then p := dxU, where u(t,x) := f~"^(i:,2;) = infj^; G [—l/2,l/2]/v(t, z) > x}. We 
have p G V2{^), and from the latter energy estimate, we deduce 

W{p) = W{p'''')- f [ \a{u{t,x))\^p{t,dx). (5.15) 
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Moreover, we have from the equation dtv{t, z) + a{z) = that the generahzed inverse u solves 
in the weak sense dfU + a{u)dxU = 0. Therefore, in the sense of distributions, we have 

dtp + d^{a{u)p) = 0. 

From Theorem 8.3.1 of [2], we deduce that p G AC f^^{{0, +oo);V2{^))- We conclude, using 
moreover flS.lSp that p is a curve of maximal slope for the functional W defined in flS.lOp . 
• Gradient flow solutions. 

This last step is a direct consequence of Theorem 11.1.3 of [2] since all assumptions of the 
Theorem have been verified above. Thus curves of maximal slope are gradient flow solutions. 
Uniqueness is obtained thanks to a contraction estimate based on a Gronwall argument as in 
|14) . This concludes the proof of Theorem 15.51 □ 

Finally we highlight that for the potential W{x) = — the problem has a special structure 
explaining why we are able in that particular case to apply the theory of gradient flows. In fact, 
in this case we have W"{x) = —6o{x) that is = in (13. 8p . Then denoting by n = W * p and 
taking the convolution by W of fl3.1ip . we deduce that u satisfies in the distributional sense 

dtu — J + w*J = 0, J = —dxA{u) + a{u)w * p. 

Therefore, in the particular case when u; = 0, this latter equation reduces to the scalar conser- 
vation law 

dtu + d:,A{u) = 0. 

Obviously, since we recover the absolute continuity of the duality solution in Proposition 15.41 
we conclude that when W = —\\x\ both concepts of solutions coincide. 
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